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Abstract. - The model of an open Fermi-system is used for studying the interplay of intrinsic chaos 
and irreversible decay into open continuum channels. Two versions of the model are characterized by 
one-body chaos coming from disorder or by many-body chaos due to the inter-particle interactions. 
The continuum coupling is described by the effective non-Hermitian Hamiltonian. Our main interest 
is in specific correlations of cross sections for various channels in dependence on the coupling strength 
and degree of internal chaos. The results are generic and refer to common features of various 
mesoscopic objects including conductance fluctuations and resonance nuclear reactions. 



Introduction. — The problem of quantum transport is generic for all realistic quantum 
systems interacting vifith environment. A transmission of a signal through a many-body 
quantum aggregate of interacting particles is essentially the main instrument in studying 
such systems and using them for practical communication purposes. Currently this is one 
of the crucial lines of development of mesoscopic physics [1] with broad applications to 
quantum information, electronics and material science. 

Historically, many ideas nowadays defining mesoscopic physics emerged in nuclear theory 
starting with Bohr's concept of reactions proceeding through compound nucleus. Low-energy 
neutron resonances in heavy nuclei present a typical example of exceedingly complex qua- 
sistationary states in an open many-body system which serve as intermediaries in reaction 
processes. Later these states provided the statistical justification for the ideas of quantum 
chaos based on random matrix theory [2] . The detailed reviews of progressing knowledge 
on quantum chaos in complex atoms and nuclei can be found in [3-6]; general features of 
mesoscopic systems of interacting fermions were stressed in [7]. The theoretical concepts of 
many-body quantum chaos are convincingly supported by the large-scale diagonalization of 
Hamiltonian matrices. Although in mesoscopic condensed matter systems one-body chaos 
often plays the main role, the interaction effects are also important and the analysis has 
many features parallel to the nuclear theory [8,9]. 
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When the lifetime of quasistationary states is getting smaU and the corresponding reso- 
nances overlap, the openness of the system becomes a decisive factor. In nuclear reactions 
this regime is called Ericson fluctuations [10], where certain fluctuations and correlations of 
cross sections are predicted. An open system can be studied by the effective non-Hermitian 
Hamiltonian [11,12] that describes the intrinsic dynamics coupled to the continuum. It turns 
out [13] that transition from isolated to overlapping resonances implies the collectivisation of 
overlapped states interacting through continuum. For a small number M of open channels, 
the restructuring of the widths leads to the segregation of M short-lived states while the 
remaining states acquire narrow widths and long lifetime. This transition is similar to the 
optical super-radiance [14]. 

One of the brightest examples of quantum phenomena in open mesoscopic systems is 
given by the universal conductance fluctuations [15, 16]. There exist well pronounced sim- 
ilarities and some differences between them and nuclear Ericson fluctuations [17]. Some 
aspects of this interrelation were studied in our previous work [18]. Using the model of 
interacting fermions analogous to the nuclear continuum shell model [19], we analysed the 
behaviour of the system in function of the intrinsic interaction strength, coupling to the 
continuum, and number of open channels. Below we study in detail how many-body chaotic 
dynamics inside the system is translated into observable features of many-channel signal 
transmission. 

Intrinsic chaos. — Our model describes n interacting fermions that occupy m single- 
particle levels of energies Eg ■ The intrinsic many-body Hamiltonian can be written as 

H^ = Ho + XV, (1) 

where Hq stands for the mean-field part describing non-interacting particles (or quasi- 
particles), and V contains the two-body interaction between the particles with the variable 
strength A. The matrix H\ of size N = ml/lnl^m — n)!] is constructed in the many-particle 
basis \k) of the Slater determinants, \k) — aj^ . . . aj^ |0), where and are the creation 
and annihilation operators. 

Each many-body matrix element Vik — {l\V\k) is a sum of a number of two-body matrix 
elements Vs-^s2S3Si involving at most four single-particle states \s). For this reason, many 
matrix elements Vik vanish, and the matrix V is very different from random matrices of 
the Gaussian orthogonal ensemble (GOE); for details, see, for example, [2,7]. The ordered 
single-particle energies, e^, are assumed to have a Poissonian distribution of spacings, with 
the mean level density l/d\, implying regular dynamics of the non-interacting system. The 
interaction V belongs to an ensemble that is characterised by the variance of the normally 
distributed two-body random matrix elements, (V^^^^^ ^^), and normalised in such a way 
that (Vil) = 1. 

It is known [7-9] that chaotic properties of the two-body random interaction (TBRI) 
Hamiltonians of the type ^ are determined by the control parameter A = X/df where df 
is the mean energy spacing between many-body states directly coupled by the two-body 
interaction. Note that d/ 3> -Do where Dq is the mean level spacing between many-body 
states. The properties of the spectra, eigenstates and observables in the model (H]) have 
been thoroughly studied in [7]. In particular, it was shown that the critical value for the 
onset of strong many-body chaos is determined by the condition A > Xcr ~ 2{m — n)/Nf 
where Nf = n{m — n) -\- n{n — l)(m — n){m — n — l)/4. 

To compare with the above model of many-body chaos, we also consider the standard 
random matrix model typically used to describe the onset of one-body chaos. The corre- 
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sponding Hamiltonian has the form, 

H^ = H°+J1Hgoe. (3) 

Here H° is a diagonal matrix with the Poissonian distribution of spacings between its ordered 
eigenvalues, and Hgoe is a. N x N matrix belonging to the GOE. In such a description, 
the Hamiltonian i?^ can be treated as a generic model describing an electron in a quantum 
dot, or as a model of optical or electromagnetic waves in a closed cavity with bulk disorder. 
The control parameter, ^ = 'jl/d^\/N , determining degree of one-body chaos is the ratio of 
the variance of matrix elements of Hgoe to the mean energy level spacing between the 
eigenstates of H° . The transition to strong chaos occurs for fi > ficr ~ 1- 

Coupling to continuum. — Our aim is to study the statistical properties of open 
systems with internal dynamics described by above two Hamiltonians. According to the well 
developed formalism [12,20,21], scattering properties of an open system can be formulated 
with the effective non-Hermitian Hamiltonian 7i, 



M 



H = H--W- W,,^\^A\A], (4) 



c=l 



where H is either iJ> or Here we neglect an additional Hermitian term (the principal 
value of the dispersion integral) that appears in the elimination of the continuum [19,22]. We 
consider the middle of the energy spectrum, where this term vanishes. The non-Hermitian 
part, W , describes the coupling between N intrinsic states Ji), Jj), through M open decay 
channels labelled as a,b,c.... The factorised structure of W is dictated by the unitarity of 
the scattering matrix. We restrict ourselves by time-invariant systems, thus the transition 
amplitudes A'^ between intrinsic states \i) and channels c are real. 

The amplitudes A'^ are assumed to be random independent Gaussian variables with zero 
mean and variance 

{A^A^) = 6,,S,,, (5) 

This is compatible with the GOE or TBRI models where generic intrinsic states coupled 
to continuum have a very complicated structure, while the decay probes specific simple 
components of these states related to a finite number of open channels (see discussion in [18]). 
Below we neglect a possible energy dependence of the amplitudes that is important near 
thresholds and is taken into account in realistic shell model calculations [19]. The effective 
parameter determining the strength of the continuum coupling can be written as 

and we consider M equiprobable channels, = 7, k'^ = k. 

All scattering properties of the system with the non-Hermitian Hamiltonian ^ are 
determined by the scattering matrix, S""^ = 6ab — i'T°'^ , with 

r'"'(i?) = f;^^('-^) 4. (7) 



h3 



The complex eigenvalues £ oi Ti coincide with the poles of the 5'-matrix and, for small 7, 
determine energies and widths of isolated resonances. In the critical region k ~ 1 with 
crossover to overlapping resonances, the width distribution displays sharp segregation of 
broad short-lived (superradiant) states and very narrow long-lived (trapped) states [13,23]. 
Correspondingly, the distribution of poles of the scattering matrix undergoes a transition 
from one to two "clouds" of poles in the complex plane of resonance energies [24] ; the number 
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of broad states coincides with the number M of open channels (the rank of the matrix W) . 
For the model ([T]), the statistical properties of resonances as a function of the interaction 
between particles and the coupling to continuum have been thoroughly studied in earlier 
papers [18,25]. 

Our main interest is in the dependence of fluctuation and correlation properties of scat- 
tering on the degree of internal chaos and strength of continuum coupling. The average 
values of reaction cross sections, a'^^{E), are fully defined by the transition amplitudes 

r°''(i;), 

a^^E) = \T''\E)\^. (8) 

In our notations the cross sections are dimensionless since we omit the common factor 
7r/fc^. We will use the terminology borrowed from nuclear physics referring to the b = a 
process as "elastic" and 6 7^ a as "inelastic" although all reactions are considered within 
the fixed energy window; in what follows we study both types of cross sections. We ignore 
the smooth potential phases irrelevant for our purposes leaving only properties due to the 
compound resonances and therefore to the intrinsic dynamics. According to the theory of 
Ericson fluctuations [10], the scattering amplitude of any process can be written as the sum 
of the average and fluctuating parts, T''''{E) = {T^^iE)) +T^''{E), with {T^^{E)) = 0. 
The average cross section, a = \T\'^ , can be also divided into two contributions, {a) = 
(cdir) + (cfl)- Here the direct reaction cross section, (adir), is determined by the average 
scattering amplitude only, while (crfl) is the fluctuational part also known as the compound 
nucleus cross section. 

Cross section correlations. — The fluctuations of both elastic and inelastic cross 
sections strongly depend on the coupling to the continuum. According to the standard 
Ericson theory, in the region of strongly overlapping resonances, k ~ 1, the variance of 
both elastic and inelastic cross sections for large M ^ 1 can be expressed via the average 
cross sections, Va.r{a) = (cta)^. Our data for the many-body Hamiltonian ^ conflrm this 
expectation for M > 10 for inelastic cross sections and any strength of interaction between 
particles. As for the elastic cross sections, a slight dependence on an internal chaos has been 
found and explained in [18]. 

Of special interest is the problem of correlations between different cross sections. The 
commonly used quantity that is discussed in nuclear and solid state physics, is the covariance 
Cfi of fluctuational cross sections, 

Cfi = (9) 

In our case without direct processes we have (cr) = (ffl), therefore, below we omit the 
subscript "fl" . The analysis of the covariance C shows that its value strongly depends on 
the type of correlations. Specifically, there are 5 types of correlations: EE - elastic-elastic, 
when both cross sections are elastic, EIi and EIq - elastic-inelastic correlations with one 
and no common channels in the scattering, and Hi and IIq - inelastic-inelastic correlations 
with one and no common channels, respectively. 

One should stress that the theoretical analysis of the covariance © encounters serious 
problems even for the GOE case in place of H in Eq. ([4]) (see also Ref. [26] for different 
model). The second term in Eq. ([9]) is defined by the second moments of the scattering 
matrix. The corresponding expressions were obtained in Ref. [21] with the use of the super- 
symmetry method. In the limit M ^ 1, they are reduced to the Hauser-Feshbach formula, 
see in Ref. [2]. It is much more difficult to evaluate the first term that is determined by 
the four-point correlation function for matrix elements. The only analytical expressions for 
this term can be found in Ref. [27]. However, the result obtained there is inconsistent with 
our numerical data, as well as with the analysis of the universal conductance fluctuations 
performed in Ref. [28]. 

In order to understand how the cross section correlations depend on the strength of 
coupling to continuum and degree of internal chaos, we performed a detailed numerical study 
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of the correlations ([9]) for two models ([T]) and All data are obtained with averaging over 
energy E at the band centre, —0.4 < E < 0.4, and over a large number of different cross 
sections belonging to one of the five groups defined above. 

In Fig. 1 we show the correlations (EE) between two different elastic cross sections, 
ct"" and a^^ with a 7^ 6, as a function of the coupling parameter. We consider here two 
limiting cases, A,/i — 0.2; 2.8, for weak and strong internal chaos, respectively. A noticeable 
dependence on the degree of chaos is clearly seen for both models. There is an excellent 
agreement between A = 2.8 and /i = 2.8 for all values of k. However, for weak coupling 
there is a small difference between the A- and /i-models. Still, the general trend of the EE 
correlations is the same for both cases. It is interesting that the symmetry between weak 
and strong coupling that is known for the GOE models, and clearly seen here for strong 
chaos, is destroyed for weak chaos, the effect lacks an analytical explanation. 
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Fig. 1: (colour online)Dependence of EE correlations on the coupling strength, k. Two limiting 
cases are shown for weak (circles) and strong (triangles) chaos. Lines and open symbols refer to the 
many-body model and full symbols to the one-body model, for = 924 and M — 4 channels. The 
average was done over 1000 realisations of random matrices and the error bars (not shown) are of 
order of the symbol size. 

Next, we consider the correlations EIq between different elastic-inelastic fluctuating cross 
sections with no common channel index, and the IIq correlations between two different 
inelastic-inelastic cross sections. These are shown in Fig. 2 for the same limiting cases as 
above, A,/Lt = 0.2;2.8. The EIq and //q correlations have a difference by a factor close to 2 
in amplitude but follow a similar trend for all three cases of internal chaos. Comparing the 
results for the two models we arrive at a similar result as before: an excellent agreement 
between A = 2.8 and fi = 2.8, and small deviations between the cases A < 1 and fi < I 
for weak coupling. This is more evident for the EIq than for the IIq correlations. The 
analogous case, when there is one common channel index between the fluctuating cross 
sections involved in the correlation function, is shown in Fig. 3. Here the EIi and Hi 
correlations differ by a factor close to 5 in amplitude but behave similarly for both limiting 
cases of internal chaos, X, — 0.2 and 2.8. Furthermore, EIi correlations are independent 
of the degree of internal chaos whereas Hi correlations become slightly smaller at stronger 
chaos. For the EIi correlations the correspondence between the two models is excellent for 
all values of k. We have to stress that the correlations for strong chaos with \, fi = 2.8, are 
in a good agreement if the Hermitian part of the Hamiltonian ([4]) is taken from the GOE. 

From our data one can see that the strongest correlations occur at perfect coupling to 
continuum, k « 1. Another remarkable property is that the correlations are either negative 
or positive depending on whether there is a common channel in the correlating cross sections 
or such channels are absent. Indeed, both EIi and Hi correlations are negative around 
K = 1, whereas the EIq and /Jq correlations are positive. 
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Fig. 3: (colour online) Dependence of the EIi and J/i correlations on k for the parameters of Fig. 
1. 



For a large number of channels, A/ 3> 1, these correlations are very weak, and they 
are ignored in the standard Ericson theory. However, they turn out to be very important 
when considering the properties of conductance fluctuations, see below. Also, in nuclear 
physics there are situations when the number of open channels is relatively small, and one 
can expect that the effect of different signs of correlations can be observed experimentally. 
One of the new applications of such experiments can be the calibration of internal chaos 
with the use of scattering data. 

Therefore, it is important to know the dependence of the correlations ^ on the number 
of channels. Our analysis has revealed that the value of the covariance C is inversely 
proportional to the number Af of terms in each group specified by the type of correlations, 
over which the averaging is performed in Eq. ([9]), see Table. Our extensive numerical data 
has confirmed this dependence, C = X/Af with some constants X that we extracted by 
fitting the data with the above dependence. This rule works perfectly starting from M = 2 
or 3. 
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Conductance fluctuations. — One of the most intriguing effects of mesoscopic physics| 
is the universality of conductance fluctuations. In order to study this effect in the framework 
of our models ^ and one should treat M/2 6-channels as left channels corresponding to 
incoming electron waves, and other AI/2 a-channels, as right channels for outgoing waves. 
Then, we define the Landauer conductance G in the standard way [15] omitting the common 
factor of 2e^/h), 

M/2 M 



(7''^ (10) 



b=l a>M/2 

The properties of the conductance are entirely determined by the inelastic cross-sections, 
b ^ a, and for equivalent channels the average conductance (jlOp reads, 

(11) 

Here T = 4k;/(1 -|- k)^ is the transmission coefficient, and F — {a^°') / {a^'') is the elastic 
enhancement factor [29] . Our data show that the value of F changes from F = 2 to F « 3 
when decreasing the strength of chaos from A,/x = oo to A,/x = 0. The last expression in 
Eq. (fTTj) is written for A/ 3> 1. As one can see, the influence of internal chaos is due to 
the enhancement factor F only. Since there is no theory relating the enhancement factor 
to the degree of chaos, we used this factor as the fitting parameter. Our data manifest an 
excellent agreement with the expression ([TT|) for various values of parameters A and /i, as 
well as M . Note that for a very large number of channels the influence of internal dynamics 
on fluctuations disappears. 

As for the variance Var(G'), the analytical results are available for the GOE case only, 
corresponding to very large values of A and ^. According to different approaches (see, 
for example, in Ref. [5]), for perfect coupling, k = 1, and very large number of channels, 
M 3> 1, this variance takes the famous values 2/15 and 1/8 for diffusive and ballistic 
transport, respectively. This result is commonly considered as a striking effect of universal 
conductance fluctuations. Our data reported in Fig. 4 clearly manifest that, in both models 
([T|) and ([2]), for K = 1 and strong internal chaos the value of Var(G') is close to 1/8. A 
small, however, clear difference from 1/8 is explained by the correction due to a finite value 
M = 20. A more general expression for a finite number of channels (for k — 1 and the GOE) 
can be found in Refs. [1, 15], 

V.riG) = 2^-M/&Ml31±2^, (12) 

and our data perfectly agree with this result. 

One can see from Fig. 4 that the variance Var(G') increases when A and ^ decrease, and 
crosses the ballistic value 2/15 close to the critical value at which the transition from weak 
to strong chaos occurs in closed models. 

One of the most interesting and new results of our study is how the internal chaos 
influences the variance of the conductance, Var(G). As one can see from Fig. 5, in the region 
with small or moderate continuum coupling, k, « 0.1 — 0.5, the value of Var(G') strongly 
depends on the strength of inter-particle interaction, A, in the model ([l]) of many-body chaos, 
and on the perturbation parameter ^ in the model ^ of one-body chaos. The strongest 
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Fig. 4: (colour online) Var(G) versus A for n = 7 , m = 14, N = 3432 and M = 20 channels for the 
many-body model HJ. 



influence of chaos occurs for k « 0.2, and the results are practically the same for both 
models, provided the appropriate normalisation of A and fi is made. This important result 
may find various applications in theory of conductance fluctuations. In particular, one may 
try to extract information about internal dynamics from experimental data when changing 
the degree of coupling to continuum (for example, the degree of openness of quantum dots) . 




Fig. 5: (colour online) Variance of the conductance for the many-body model ([3}, left panel, as 
compared to the one-body model, right panel. The average was done over 700 realisations of 
random matrices with = 924 and M — 10. 



An important feature of the data in Fig. 5 is that for k « 1 the influence of intrinsic chaos 
is weak. This can be explained as follows. When coupling is perfect, an interaction with the 
continuum through forming broad states in various channels is very strong in comparison 
with an internal process of chaotization, therefore, the latter may be neglected. We found 
that for K 1 the sensitivity of the conductance fluctuations to the degree of chaos decreases 
with an increase of number of channels M . 

It is instructive to show that the main properties of conductance fluctuations cannot be 
explained if one neglects cross section correlations discussed above. For the first time, the 
role of these correlations in application to conductance fluctuations has been discussed in 
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Refs. [30-32]. The variance Var(G) can be rewritten as 



Var(G) = Ml ( 



T 



) 



2 



F + M -I 



(13) 



where A''! = L(M-2), A^o = L{L-M+l), L = M2/4 and the terms Ci and Co stand for the 
III and IIq correlation functions, respectively, see above. If one neglects the correlations, 
the first term gives 1/4 (for T = 1 and Af 3> 1), instead of 1/8. Our analysis shows that for 
the strong interaction and M ^ 1 one obtains Ci « —M~^ and Co ~ 2A/~^. Therefore, in 
the limit of a large number of channels, the first term in r.h.s. of Eq. (|13p that equals 1/4, 
is cancelled by the second term, and the third term tends to 1/8 resulting in Var(G) = 1/8. 

This result clearly demonstrates the crucial role of correlations determining the conduc- 
tance fluctuations (see, also, Ref. [28]). Remarkably, the //i-correlations cancel the first 
term 1/4, and the value 1/8 is due to the //o-correlations (term Cq) only. A highly non- 
trivial role of the Ci and Cq terms can be also manifested in the correlations of speckle 
pictures [31]. 

Conclusions. — To conclude, we studied the interplay of complicated intrinsic dy- 
namics and coupling to the outside world for typical quantum systems with one-body or 
many-body intrinsic chaos, the first one coming from the disordered single-particle spec- 
trum and the second one emerging as a result of inter-particle interactions. The openness 
of the system is described by the effective non-Hermitian Hamiltonian that fully respects 
the unitarity requirements and allows to calculate, in the same framework, cross sections of 
various processes, their fluctuations and correlations. As a manifestation of the general fea- 
tures of underlying physics, the models are equally valid for description of nuclear reactions 
with the transition from isolated to overlapping resonances and for conductance fluctuations 
in mesoscopic condensed matter devices. We found that the correlations of inelastic cross 
sections are very different for the processes with and without a common channel, being 
negative in the latter case in the region of perfect coupling when the typical decay widths 
and resonance spacings are of the same magnitude. Another important result is that for 
the conductance fluctuations the dependence on the degree of intrinsic chaos is strong at 
intermediate continuum coupling, in contrast with the region of perfect coupling, for which 
the continuum dominates, and the fluctuations are known to be independent of internal dy- 
namics. Many results of the conductance theory are numerically confirmed being explained 
by the specific correlations of partial cross sections of very general origin. 
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